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Abstract 

The extension of the statistical parton distributions to include their trans- 
verse momentum dependence (TMD) is revisited by considering that the pro- 
ton target has a finite longitudinal momentum. The TMD will be generated 
by means of a transverse energy sum rule. The new results are mainly rel- 
evant for electron-proton inelastic collisions in the low Q 2 region. We take 
into account the effects of the Melosh-Wigner rotation for the helicity distri- 
butions. 
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1 Introduction 



In 2002 [TJ we have proposed a description of the parton distributions inspired 
by quantum statistics with some robust phenomenological motivations. 

1) The defect in the Gottfried sum rule [21 E], implies d(x) > u(x), ad- 
vocated many years ago as a consequence of the Pauli principle [I]. This 
inequality has been confirmed by Drell-Yan production of muon pairs [5], up 
to some moderate x values and remains to be verified for x > 0.2 or so. 

2) The dramatic decrease at high x of the ratio of the structure functions 
F^x) / F^x) [6], strongly connected to the behavior of d(x)/u(x), which is 
predicted to flatten out for x > 0.6 [7]. 

3) The fast increase at high x of the double longitudinal-spin asymmetry 
A±(x), which implies the dominance, in that region, of the w-quark with the 
helicity along the polarization of the proton target. 

The correlation between the shapes and the first moments dictated by the 
Pauli principle, allowed us to describe with a small number of parameters 
both the unpolarized and polarized distributions with the important conse- 
quence to predict a positive value for A«(i) and a negative one for Ad(x), 
giving a positive contribution to the Bjorken sum rule [8]. 
Typical predictions of our approach, as the monotonic increase of the posi- 
tive ratio Au(x)/u(x) and decrease of the negative ratio Ad(x)/d(x), are in a 
good agreement with experiment [HI ED]- It should be stressed that our pre- 
diction on Ad(x)/d(x) does not correspond to the belief that this quantity, 
according to the counting rules, should change sign to reach +1 for x — 1. 

In earlier works [TT| IT2] . we succeeded to explain the arbitrary factors, 
which were necessary to agree with data by the extension to the transverse 
degrees of freedom and we constructed a set of transverse momentum depen- 
dent (TMD) statistical parton distributions. Here we want to improve our 
approach by considering the longitudinal momentum of the proton target P z 
finite and not oo, as it is only in the limit Q 2 — > oo. 

The paper is organized as follows. In the next section we review the construc- 
tion of the statistical distributions and then we derive the exact transverse 
energy sum rule, the TMD statistical distributions and the expression of P z as 
a function of Q 2 and x. In section 3 we recall the expression of the Melosh- 
Wigner transformation which has an important role for the TMD helicity 
distributions, as will be shown in section 4. We give our concluding remarks 
in section 5. A useful transformation to simplify the calculation of integrals 
over the transverse momentum of the partons is given in the Appendix. 
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2 The statistical approach 



In our previous papers [TJ [7] we took the following expression for the non- 
diffractive contribution of the quark distributions q h (x) of flavor q and helicity 
h, at an input energy scale Ql 

exp[(a; - X^ q )/x\ + 1 

where X$ q plays the role of a thermodynamical potential and x of the universal 
temperature. Correspondingly for antiquarks we took the expression 

Ax^ 

^ = X^[eM(x + X 0q h )/x] + l] ' (2) 

which shows a strong connection between quarks and antiquarks of opposite 
helicity. In order to take into account the uniform rapid rise of all distri- 
butions in the very low x region, we had to add to the above nondiffractive 
contributions, for quarks and antiquarks, the following diffractive contribu- 
tion 

Ax^ 

XqD{x) = eMx/x) + V (3) 

which is flavor and helicity independent. 

Finally the expression for the gluon distribution was 

xG(x) = , (4) 

exp{x/x) — 1 

since the potential of the gluon Xqg must be zero. 

From the well-established features of the u and d quark distributions, ex- 
tracted from Deep Inelastic Scattering (DIS) data, we could anticipate some 
simple relations between the corresponding potentials, namely 

(5) 
(6) 

+ X o~ d • (7) 

Clearly these partons distributions pi = q,q,G must obey the momentum 
sum rule, which reads 

/ xpi{x)dx = 1 . (8) 
i Jo 



X 0u 


> 


X 0u 


X 0d 


> 


X 0d 


X 0u + x ou 


> 


X 0d 
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To account for the factors Xq in the numerator of Eq.([T]) and X^f 1 in the 
denominator of Eq.([2]), we extended our analysis to the transverse degrees 
of freedom [11] . If pt stands for the parton transverse momemtum, the 
corresponding TMD statistical distribution pi(x,p^) is normalized such that 
J pi{x,p\)dp\ = Pi(x). Let us denote by E the energy of the proton target of 
mass M and longitudinal momentum P z . By definition the z axis is along the 
direction of the proton momentum and this implies its transverse momentum 
Pt to be zero. The deep inelastic regime, where parton model may be applied, 
is at high Q 2 and in that limit, one may neglect M with respect to P z . 
However at very small x, one may not neglect pt with respect to xP z . The 
proton energy is taken to be E = P z + M 2 /2P Z . Similarly the energy of a 
massless parton, such that xP z » pt, is xP z + p T /2xP z . So by writing a 
sum rule for the energy of the partons analogous to ([8]), we see that it does 
not depend on P z and one gets now 



As explained first in Ref.[TT] and more recently in Ref.|12j. an improved 
version, this energy sum rule implies the p\ dependence of the nondiffractive 
contribution of the quark distributions to be 



where Yq is the thermodynamical potential associated to the parton trans- 
verse momentum px and l//i 2 is a Lagrange multiplier, whose value is deter- 
mined by the transverse energy sum rule. 

Before moving on, we would like to make an important remark on the im- 
plication of the above transverse energy sum rule. If one makes the com- 
monly used assumption that Pi(x,p T ) obeys a factorization property, such 
as, pi(x,p T ) = Pi(x) ■ fi{p T ) } then Eq.(j9]) implies that dx pi(x)/x must 
converge. Consequently for all partons, the small x behavior is strongly 
constrained, since it must vanish like x a \ with a« > 0. This is in obvious 
disagreement with the observed rapid rise of sea quarks and gluon in the very 
small x region. We conclude that our transverse energy sum rule does not 
allow the simplifying factorization assumption for the TMD distributions. 
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(10) 



exp(p T /x[i 2 - Yq ) + 1 ' 
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2.1 Transverse energy sum rule and px dependence 

The approximations P z » M and xP z » pr, valid in the limit P z = oo, 
bring in an unphysical singularity for x — > 0, where evidently xP z is not 
larger than px- We have the relations 

E = J> (11) 

E-P z = ^2(Ei- Pgi ), (12) 

i 

where p 2 j = £jP 2 , £7j = a/ x 2 P 2 + p 2 ^ + mf is the energy of quark i, and 
E = y/P* + M 2 . From now on we shall neglect the quark mass r?v By 
multiplying both sides of Eq. (ll2j) by E + P z and inside the summation, both 
numerator and denominator, by E^ + p Z i, we get 

M 2 = (E + P.) ]P {Ei ~^[ Ei + Pzi) p^pl,) (13) 
*Y {Ei + p zi ) 

= ^ + ^E7t A p ttO*,&)- (14) 
i V ^ p 2 + Pt; + XiPz 



Therefore instead of Eq.([9]) if we consider Eq.l jTjj) in the continuum limit, 
the exact transverse energy sum rule reads 

M2 = (VpFTW + P .) E / ix [ ^ - - (15) 



since P z >> M. 



By comparing the integrants of Eq.Q and of Eq. ffl6l) . we see that the 
Pt dependence of the parton distribution, which was driven in Eq. ffTO]) by 
Pt/x/j, 2 , is now replaced by the following expression involving P z 

' (17) 



exp[2p 2 T /[^ 2 (x + ^x 2 +p 2 T /P 2 )} - Y h q ] + 1 ' 
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As a result the TMD of the nondiffractive contribution of the quark distri- 
butions are now given by the following expression 



xq h (x,PT, P z ) — r, vh\ /-i , 1 x 
exp[(:r - X(} q )/x] + 1 



exp[2p£ /[^{x + y/x 2 +p 2 T /P 2 )} - Y Q \\ + 1 ' 

where F q is a normalization function , which has the dimension of the inverse 
of an energy square, and which will be explained in more details below (see 
Section 4). 

2.2 Expression for P z finite 

The value of P z for finite Q 2 , which seems more natural to us, is its value 
in the reference frame where the final hadrons are at rest. In this reference 
frame the proton target and the virtual photon have opposite momentum. 
To go to this frame from the lab frame where the proton is at rest, we have 
therefore to consider the Lorentz transformation with velocity 

^ = ^-T7> (19) 
where (go, ~<t) is the virtual photon momentum. In that reference frame 

(20) 



M 3* M\?\ 



l^l 2 y / q 2 + 2Mq + M 2 - [f\ 



which implies 



M 2 {Q 2 + q 2 ) 

2 2Mq + q 2 + M 2 - \^\ 2 ' 1 1 

Now from the definition of the scaling variable x = Q 2 /2P ■ q, we get 

Q 2 = 2xMq , (22) 

so we can write 



m 2 (Q 2 + Jm.) _ xM 2 (J w + i) 

31 + M 2 -Q 2 l-x + ^- 



(23) 
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Finally we get the following simple expression of P z in terms of x and Q 2 




10 3 10 3 10 1 1 

X 

Figure 1: P z versus x for different Q 2 values. 

notice that when Q 2 = oo, then P z = oo if x ^ 0. However for Q 2 finite, 
when x — > 0, we also have P z — > oo. 

From Eq. f )24p one realizes that the effect of the finite value for P z in Eq. (j!8p 
is relevant at high values of p?n, for x in the neighborough of 1 and for small 
values of Q 2 . 



3 The Melosh-Wigner transformation 

So far in all our quark or antiquark TMD distributions, the label " £ /i" ' stands 
for the helicity along the longitudinal momentum and not along the direction 
of the momentum, as normally defined for a genuine helicity. The basic effect 
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of a transverse momentum 7^ is the Melosh-Wigner rotation [T3"| 11%] , 
which mixes the components g 1 * 1 in the following way 

q +MW = cog 2 q q+ + gin 2 g g - g -MW = cog 2 g ? - + gin 2 £ ?+) ( 25 ) 

where, for massless partons, 

6 = arctan ( — — — ) , (26) 

PO+Pz 



with p = ^/pt + pI- 

Consequently q = q + + q~ remains unchanged since q MW = q, whereas 
we have 

Aq MW = (cos 2 6 - sm 2 6)Aq . (27) 

The angle of the Melosh-Wigner transformation Eq.( )26|) . originally derived 
in Ref . [13] , has the correct property to vanish when either px = or 
P z -»■ 00. 

From simple calculations @ we get 

2n ■ 2n (Po + Pzf ~ Pt Pz xP z 

cos — sin 



(Po + Pz) 2 +Pt Po ^x 2 P 2 + p\ 
1 1 1 



53^ V ^? " 

where we have used the variables £ and 77 introduced in the Appendix and 
Eq.(0O]). 



In the last term of Eq. (l28j) by using the expression of £ given in Eq. flA.7j) 

2 2 

we obtain 1/(1 + -^ps) and this term cancels 1 + in Eq. flA.9[) . 



The consequences of this cancellation will be discussed in the next section. 

4 Some explicit expressions for the TMD dis- 
tributions 

One of the goals of this section is to show how to use the normalization 
function introduced in Eq. flTBl) and then how to evaluate it. 



1 The angle we used in Ref. [12], which was phenomenological, didn't have this property. 
2 This result has been also used in several papers on the proton spin puzzle |15[ 116] . 
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Following Eq. OA.ip and the results of the Appendix, the nondiffractive 
contribution of the quark distributions q h (x) read 



xq\x, p 2 T , P z )d V \ = pthJ^Ll f x ln (1 + e YZ q)[l + Rh] 
o expl(x-X$ q )/x] + 1 

(29) 

For illustration we display in Fig. 2, using the parameters /i 2 = 0.198GeV 2 
and = 1.122, determined in Ref. [12J. The other correction factors R~, 
R^, are very similar, so it is clear that the effect of these correction fac- 
tors is limited to the small Q 2 region. By using Eq. (JTJ) one gets the following 



Q = 2GeV 



10" 



10 



10 



; Q = lOGoV 



; Q = lOOGeV" 



10" 



10" 



10 



10 



Figure 2: i?+ versus x for different Q 2 values. 



condition for F q 



F q h (x,P z ^ 2 ) 



AX* q x b 



fi 2 xln{l + e Y ^)[l + R^] 



(30) 



The parameters A, b, X^ q are known, P z is a function of Q 2 , x, the only 
parameters to be fixed are fi 2 and Yq q . It is worth noting that since q h (x) 
should not depend on P z , F q is such that it will absorb the P z dependence 
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occuring only in the correction factor Bh. 

The unpolarized TMD distributions are defined as 



xq(x,PT,P z ) 



ex P [(x - X+)/x] + 1 exp[2j%/\p?(x + ^x 2 + p* T /P?)]-Y+] + l 

1 F-(x,P z ,fx 2 ) 

exp[(x - X 0q )/x] + 1 exp[2^/[/i 2 (x + y/x* + p 2 T /P 2 )] -Y c 



(31) 



OqrJ 



For illustration we display in Fig. 3, the TMD of the nondiflr active con- 
tributions of the quark distributions xu(x,p^,Q 2 ) and xd(x, p?, Q 2 ), using 
the above normalization condition Eq. (!30|) and the parameters of Ref.|12j. 
It is clear that all these pr distributions are close to a Gaussian behavior, 
but with a x-dependent width. As expected, the fall off in p T is faster for 
smaller x values and we see that the effect of the correction factor R^, which 
decreases the distributions, is more important near p? = 0. 




0.2 0.3 
Pt (GeV) 



0.2 0.3 
Pt (GeV) 



0.5 



Figure 3: The TMD nondiffractive contribution of the quark distributions 
versus pr, at Q 2 = 2GeV 2 , for two x values: solid and dotted lines x = 0.01, 
dashed and dot-dashed lines x = 0.1. (upper curves Rh = 0, lower curves 
R\ ^0) 
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Let us now turn to the diffractive contribution Eq.Q. Since b < 0, one 
cannot introduce the px dependence similarly to the nondiffractive contri- 
butions, because it generates a singular behavior in the energy sum rule, 
when x — > 0. In order to avoid this difficulty, as in Ref. \T'2\, we modify our 
prescription by taking 

xq D (x,p 2 T ,P z ) = ) x 

exp[x/x\ + 1 

F D (x,P„y?) 

exp[2p 2 T /[^x(x + ^x 2 +p 2 T /P 2 )]] + 1 ' 

whose pr fall off is stronger, because /i 2 is now replaced by xfi 2 . In order to 
make a link with Eq.(j3]), we have to compute the integral 

✓•oo 

xq D (x)= / xq D (x,p 2 T ,P z )dp 2 T . (33) 
Jo 

By using Eq. (lA.lOp with the substitutions Yq = and /i 2 — > xfi 2 , we easily 
deduce from Eq. f[3"31 

Ax b ~ 2 

since Li 2 (— 1) = — 7r 2 /12. This result is similar to what was found in Ref.|12j. 
because the correction factor is very small. 

Concerning the TMD gluon distribution, it reads similarly 

xG(x, Pt, P z ) = : t— : - 

exp|x/x| — 1 

F G (xP„ii 2 ) 

ex V[ 2p 2 T /[fi 2 x(x + ^Jx 2 + p 2 T /P 2 )] + Yq] — 1 

Here we must introduce a negative potential —Yq to avoid a singularity when 
Yq — > 0. In order to make a link with Eq.Q we have to compute the integral 



which reduces to 

"OO J 2 



xG(x) = I xG(x,Pt, P z )dp% , (36) 
dp\ 



o exp [2p 2 T /[fi 2 x{x + ^x 2 +p 2 T /P 2 )} + Yq] — 1 

- /i 2 *(l + ^)dr/ 



exp[?7 + Yq] — 1 ' 



(37) 
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and by using Eq. 

F G (x,P z ,^[ 



one obtains 



(3? 



A G x bG ~ 



ln (1 - e-yo) + ^a n 2 _ iyj + y G \n (e^ - 1) + Li 2 (l - c *b)] 



This result is similar to what was found in Ref.|12j. because the correction 
factor is very small since we take Yq = 10~ 6 . 

Before we move on, we would like to compare our results with the relativis- 
tic covariant approach [17J, where they introduce the variable x + p^/xM 2 
combining the x and dependences. We show in Fig. 4 the result of our 
calculations which can be compared with the results of Ref.[PT] displayed 
in their Fig. 1. The two results are compatible with a broader shape for 
increasing x, but one notices that the px fall off is less rapid in our case. 




p T /M 



p T /M 



Figure 4: The TMD distributions for u and d quarks, versus pr/M, at Q 2 = 
4GeV 2 , for two x values: solid lines x = 0.18, dashed lines x = 0.3. 

The TMD helicity distributions xAq(x,p^, P z ) will be defined as above, 
by substracting instead of adding the two helicity components, so there is no 
diffractive contribution. 

It is clear that after integration over p?, Aq(x) does not depend on P z . 
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Figure 5: The u and d quark helicity distributions versus x: xAq(x) (dashed 
line) and xAq MW (x) (solid line) 

However this is not the case for the helicity distributions modified by the 
Melosh-Wigner transformation, since after taking into account the cancella- 
tion obtained in Section 3, we have 

xAq MW (x,P z ) = 

—i- y + P z , fi 2 )fi 2 x In (1 + e *3) 

exp[(x - XfiJ/x] + 1 

r? 1-\ /-I -l. 1 Pz ' /i 2 ) Aln (1 + e y °«) . (39) 

exp[(z - A 0g )/xJ + 1 

By using Eq.( )30|) one finds 



[exp[(x-X +)/a;] + l][l + J R+] 
Ax 6 

"[exp[(x-X -)/x] + l][l + J R ? 



(40) 



As expected, in the limit P z — > oo, the Melosh-Wigner transformation be- 
comes an identity, so xAq MW (x, P z ) — > xAq(x). For illustration we display 
in Fig. 5, xAq(x) and xAq MW (x) for Q 2 = 2GeV 2 , which shows the effect 
of the Melosh-Wigner rotation, mainly in the low x region. 
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It is interesting to note that \Aq MW (x)\ < \Aq(x)\, as expected from 
some earlier work [To] . 

5 Concluding remarks 

We have presented the extension of the statistical parton distributions to 
include their transverse momentum dependence, by considering that the pro- 
ton target has a finite longitudinal momentum P z . This situation generates 
some correction factors, which are only relevant in the small Q 2 region and 
for rather low x values. The TMD distributions were generated by means 
of a transverse energy sum rule which implies that a simplifying factoriza- 
tion assumption is not allowed, as explained above. This sum rule has been 
used by other authors [18] who have, as in the present work, non factor- 
ized TMD distributions, as well as in the relativistic covariant approach [T7] . 
The TMD diffract ive part of the quark (antiquark) distributions q (z,Py) 
and the TMD gluon distribution G(x,pji) had to be treated in a different 
way, as in Ref.[T2], in order to avoid a singularity in the energy sum rule. As 
a result their contributions to the sum rule are very small and the nondiffrac- 
tive contributions dominate. Our approach involves the parameter /z 2 , which 
plays the role of the temperature for the transverse degrees of freedom, whose 
value is determined by the sum rule, as given in Ref.[T2] which is probably 
an upper bound. If the TMD gluon distribution contributes significantly to 
the energy sum rule, as it does to the momentum sum rule, we might obtain 
a much smaller value, but this remains to be proven and will be the subject 
of a future work. 

We have also shown the importance of the remarkable Melosh-Wigner 
transformation, whose effects are significant only for a finite P z . 

At the difference of other approaches, where statistical concepts are used 
in the target rest frame and then a boost is applied to a large P z frame, we 
made the choice to consider directly the large P z frame and the x variable to 
define and the shapes of the different distributions.. 

There is a phenomenological evidence that the partons dominating the 
large x regions are not the same, as those which dominate the low x region, 
as one could mainly find by boosting an isotropical rest frame distributions. 
So we think that there is a deep theoretical reason to settle the statistical 
concepts directly in the x variable related to the foundation of statistical 
mechanics. One may also cast some doubts on the use of the statistical 
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approach, since the total number of valence partons is small, 2-u quarks and 
1-d quark in the proton, but the fact that one writes a probability density 
makes reasonable to apply a statistical approach to probabilities. 
Finally, we would like to refer to an interesting recent paper suggesting a 
duality principle between our approach and a thermal description of the PDF 
[19] . This duality allows them to introduce an effective temperature T~120- 
150MeV, which is approximately the same for longitudinal and transverse 
momentum. The comparison with our results may be done in the Boltzmann 
limit, where we neglect the effect of quantum statistics, which is crucial to 
get the phenomenological successful shape-first moment correlation for the 
valence partons and the isospin and spin asymmetries of the sea. In that limit 
we get for the longitudinal temperature T) = Mi/2, which gives 47MeV. To 
get the transverse temperature T t , always in the Boltzmann limit, one should 
follow the method described in Ref.[20] leading to T t = fj,y/E/2, which gives 
70MeV, not too far from T\. As noted in Ref.|19j. there is no contradiction 
in getting different values in our approach and in the thermal description of 
the PDF. 



A Appendix 

Let us recall that the nondiffractive contribution of the quark distribution 
q h {x) introduced earlier, must be obtained by an integration over p?, as 
follows 

roo 

xq\x) = / xq h (x,p 2 T ,P z )dp 2 T , (A.l) 
Jo 

where xq h (x,p^, P z ) is given by Eq.( fl8|) . So if we introduce the variable 
£ = Pt/^ 2 , we have to consider the following integral 

X ^ (A.2) 



o exp[ - 7 A_ -y*] + l 



A simple inspection of Eq. flA.2j) shows that the pt dependence is a compli- 



cated expression which involves a square root, so the integral is certainly not 
tractable analytically. One way to by-pass this difficulty is to transform the 
integrant in the form of a usual Fermi-Dirac function and consequently to 
compute the associated differential element. In the above Eq. (1A.2j) let us 
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perform the change of variable 



V = i=== , (A.3) 



which can be written as 



- = A/ 1 + i ^ + l- ( A - 4 ) 
77 y xrt: 



Now from Eq.( 1A.4l) one has 



xP 2 r/ 7] 

and finally the relation between 77 and £ 



1 + ^ = ^-H + i (A.5) 

x p2 ^2 ^ 

-(- - 1) , (A.6) 



Z = V + ^. (A.7) 



AxP* 

By differentiation one gets 

d£=(l + ^jk)dv, (A.8) 

so with the above transformation, the integral (1A.2|) can be rewritten in a 
simplified form, close to a Fermi-Dirac distribution, 



l exp[r) - Y h q ] + 1 
To summarize we have shown that 



(A.9) 



exp[2 P y[^(x + ^xZ+rt/P?)} - Y h q ] + 1 
r°° xu 2 (l + -^n)dr] , ,,2 

J expfo - Y h q ] + 1 1 1 J 2xP2 21 n 

H 2 x In (l+e y &)[l+i#], (A.10) 
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where Li 2 denotes the polylogarithm function of order 2 and Eh is a correction 
factor i?g = — ^pi—^-^pP-, which involves the P z dependence. 

z ln(l+e °i) 

We see that in the limit P z — > oo, we recover the original formula given in 
Ref.[12], for the integral of the TMD, since the correction factor disappears. 
In DIS, one can neglect xM 2 with respect to Q 2 , so P 2 = Q 2 /4x(l — x) and 



r>h 2^(1 -x) Li 2 (-e y A) f , 

Rq ~ Q 2 ln(l + e^) ' (A ' n ' 



and we note that > 0, because Li 2 (— e °«) < 0. 
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